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The decay widths of the d* —¥ dn°n 0 and d* —¥ d-K + n~ processes are explicitly calculated in terms 
of our chiral quark model. By using the experimental ratios of cross sections between various decay 
channels, the partial widths of the d* —► pmr°n 0 , d* —> pnn + n~, d * —> ppn 0 TT^, and d * —> nmv + Tv° 
channels are also extracted. Further including the estimated partial width for the d* —¥ pn process, 
the total width of the d * resonance is obtained. In the first step of the practical calculation, the effect 
of the dynamical structure on the width of d* is studied in the single AA channel approximation. 
It is found that the width is reduced by few tens of MeV, in comparison with the one obtained 
by considering the effect of the kinematics only. This presents the importance of such effect from 
the dynamical structure. However, the obtained width with the single AA channel wave function 
is still too large to explain the data. It implies that the d* resonance will not consist of the AA 
structure only, and instead there should be enough room for other structure such as the hidden-color 
(CC) component. Thus, in the second step, the width of d* is further evaluated by using a wave 
function obtained in the coupled AA and CC channel calculation in the framework of the Resonating 
Group Method (RGM). It is shown that the resultant total width for d* is about 69 MeV, which is 
compatible with the experimental observation of about 75 MeV and justifies our assertion that the 
d* resonance is a hexaquark-dominated exotic state. 
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I. INTRODUCTION 

In recent years, the CELSIUS/WASA and WASA@COSY Collaborations successively reported the observation of 
a resonance-like structure in the double pionic fusion channels pn -A dTr°Tr a and pn —> dn + Tr~ when they studied 
the ABC effect and in the polarized neutron-proton scattering [ 1 H 3 . They mentioned that because the width of the 
structure is rather narrow, which is three more times smaller than 2 Ta in the conventional A A process, the observed 
data cannot be explained by the contribution from either the Roper excitation or the t-channel AA process. Therefore, 
they proposed a d* hypothesis, in which its quantum number, mass and width are I{J P ) = 0(3 + ), M r: 2.36 GeV 
and r ft: 80 MeV [I] (in their recent paper Q, they take averaged values over the results from elastic scattering and 
two-pion production, i.e. M r; 2.375 GeV and T r; 75 MeV), respectively, to accommodate the data. Because “the 
structure, containingsix valence quarks, constitutes a dibaryon, and could be either an exotic compact particle or a 
hadronic molecule” 0, this result causes physicists’ special attention. 

In fact, the existence of the non-trivial six-quark configuration with I(J P ) = 0(3 + ) (called d* lately) has intensively 
been studied since Dyson’s estimation @. A variety of methods or models, such as group theory [§], bag quark 
model [?], quark potential model |M3, etc., have been employed to investigate the structure of d*, among which 
even some investigations produced a mass close to the recent data, they are either not a dynamical calculation, or a 
calculation without the width prediction or with an incorrect width prediction. It should specially be noted that in 
one of those papers 0, one performed a coupled channel dynamical calculation in 1999 where a AA channel and a 
hidden-color channel (denoted by CC) are included and the predicted mass is about 40 — 80 MeV. This means that 
in this structure, there might exist a six-quark configuration, which coincides with COSY’s assertion. Nevertheless, 
in that paper, the width of the state has not been calculated. 

After COSY reported their finding, many investigations have been devoted to this aspect. There are mainly three 
kinds of models on the structure of the d* resonance: a) It is a AA resonance [Tl|] . The authors in Ref. El performed 
a multi-channel scattering calculation and obtained a binding energy about 71 MeV with respect to the AA threshold 
and a width about 150 MeV where Tnn — 14 MeV and Ti n ^ = 136 MeV. b) It is dominated by a “hidden-color” 
six-quark configuration. Bashkanov, Brodsky and Clement [121 ] argued in 2013 that this hidden-color structure is 
necessary for understanding the strong coupling of d* to AA. Later, Huang and his collaborators made an explicit 
dynamic calculation in the framework of the Resonating Group Method (RGM) 0] and showed a binding energy 
of about 84 MeV and almost 67% of “hidden-color” configuration in d*. This imp lies that d* is probably a 6-quark 
dominated exotic state, c) It is a result of the ANn three body interaction 0]. In order to justify which one of 
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these three is more reasonable, a detailed calculation, especially the decay width, should be performed and further 
experimental investigation should be carried out. 

In this paper, we focus on d* width study. We would firstly exam the effect of the dynamical structures of the 
d* and deuteron bound states on the decay width of d*, and consequently fetch out the contribution from the A A 
structure of d* with J p = 3 + . Then, we would estimate the total width of d* by including the contributions from 
other possible decay channels. At the beginning, we temporarily assume that d* is composed of the AA structure 
only. In the calculation, the extended chiral SU(3) quark model is employed, because this constituent quark potential 
model can successfully reproduce the spectra of baryon ground states, the binding energy of the deuteron, the nucleon- 
nucleon (NN), Kaon-nucleon (KN) scattering phase shifts, and the hyperon-nucleon (YN) cross sections (for details 
see Refs. |15l - H7| ). With the same set of model parameters fixed in explaining the above mentioned data, the bound 
state problem of the AA system is solved and the realistic wave functions of d* and deuteron are obtained via the 
dynamical RGM calculation. With these wave functions, the two-pion decay width of d*(2380) in the process of 
d* —> dirir is calculated on the quark level, where the chiral effective Lagrangian of quark-quark-pion is employed. In 
terms of the experimental data of other observed decay channels such as d* —> np7r°7r°, d* — > npn + Tr~, d* —> nmr + Tr °, 
d* —>■ pp7r°7r~, etc., the total width of d* is estimated, and the role of dynamical structures to the decay width is 
analyzed. The result with the single AA channel assumption exhibits the importance of the dynamical structure effect 
which reduces the decay width by about few tens of MeV. However, the width is still larger than the experimentally 
observed value, so that the other structure in d* should further be considered. Subsequently, we evaluate the width 
of d* with the wave function obtained in the coupled A A and CC channel RGM calculation |k|]. The resultant total 
width of d* is about 69 MeV, which is compatible with the experimental data. In the next section, the formulism is 
briefly given. The numerical results and discussion are presented in the final section. 


II. BRIEF FORMULISM 


Referring to Ref. |l8j, the phenomenological effective Hamiltonian for the quark-quark-pion interaction in the 
non-relativistic approximation is 

V = 9wt 54r^x (2jr)3W (1) 

where g qq7 r is the coupling constant, <j> stands for the n meson field, and k n are the energy and three-momentum of 
the 7r meson, respectively, and a (r) represents the spin (isospin) operator of a single quark. The wave functions are 


N>= V2 Xp ^ P 


XaV’a 


®n(p, A) 


for the nucleon and 


A > = Xs^s^a(p, A) 


( 2 ) 


(3) 


for the A(1232) resonance. In Eqs. (2-3), x and ip stand for their spin and isospin wave functions, $jv(p, A) and 
are the spatial wave functions of the nucleon and A resonance, respectively, and p and A are the Jacobi coordinates 
for the internal motion. Then, the decay width for A —► ttN reads 

(4) 


3n 


where uj n ,N = \J + k% are the energies of the pion and nucleon, respectively, ~ 0.229 GeV, and I a denotes 
the spatial overlap integral of the internal wave functions of the nucleon and the A resonance. By fitting the measured 
width of 117 MeV for A 3 / 2 +(1232) Qjj], one gets G = g qqn I 0 ~ 5.41 GeV -1 which is the product of the coupling 
constant g qqn and the spatial integral Iq. 

Now using the knowledge of AIa-m tN obtained above, we can estimate the decay width in the d* —> dnir process. 
The transition matrix element between the initial state d* and the final state d7r°7r° can be written as 


d 3 q 


Xd(Q~$ 12 ) . Xd(<l+?k 12 ) 


\-E A (q) - E N (q - £q) - uq 
x* d (q + \k 12 ) 


E A (q) - E N (q - k 2 ) ~ w 2 
Xd(q- 5 * 12 ) 


E A {~q) - E N (-q - ki) -uq E A (-q) - E N (-q - k 2 ) ~ W 2 -I 


Xd*(q), 


(5) 
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where i and / stand for the initial d* state with quantum numbers ((Sms) = (3m<j*)) and the final deuteron state with 
((Sms) = (1 m d )), respectively, is the spin (isospin) factor shown in the appendix, = F(k\ 2 ) = ^ 27r p/^= , 

fci 2 = ki — ^ 2 , Wi t 2 = + fc \ 2 . Xd(oO and Xd-(q) are, respectively, the relative wave functions of the final deuteron 

(between the two nucleons) and the initial d* (between the two-As) where q = \(p\ + P 2 + P 3 ~ P 4 — P 5 — Pe) with p) 
being the momentum of the i-th quark. Four terms in the bracket of Eq. ([5]) are related to the propagators of four 
sub-diagrams in Fig. 1. 





FIG. 1: Four possible emission ways in the decay of the d* resonance composed of the AA structure only. Two pions with 
momenta of k\^ are emitted from one of the three quarks in 2 As, respectively. 

With the transition matrix element _A/f^ 7r °, the decay width of d* in the d* —> dirir channel can be evaluated by 

J d 3 k 1 d 3 k 2 d 3 p d (2Tr)S 3 (ki + k 2 + p d )6(uj kl + uj k2 + E Pd - M d .) \ M*/ | 2 , (6) 

where uj kl) k 2 are the energies of the two outgoing pions, E Pd is the energy of the outgoing deuteron with momentum 

Pd, and the bar on the top of the transition matrix 71 means that this matrix element is averaged over the initial 
states and the summed over the final states. The factor of 2! is due to the two identical pions in the final states. 

In the practical decay width calculation, one needs the explicit deuteron and d*(2380) relative wave functions. These 
wave functions can usually be taken from the realistic solutions of the system considered. In this work, we obtain 
them by dynamically solving the RGM equation in the extended chiral SU(3) quark model fl3j where the binding 
energy of deuteron is e = 2.2 MeV and the binding energy of d* is e ~ 62 MeV in the single A A channel approximation 
and e ~ 84 MeV if the CC channel is further considered, and consequently, the mass of d* is M d * = 2Ma — e. In the 
coordinate space, the wave functions of the deuteron and d* systems can also be expressed, respectively, as 

d = [ 4 >n(^1:^2) 4>n(U^ 5) Xd(R) ] C(S/) = ( 10 ); (V 

’I'd* = [ </>a(€i,&) <M(£4,£5) XAa(R) + 0 c (£l 162) ^(^4,^5) XCc(R) ] C(S/) = ( 30 ); 

where 0 a, 0c denote the internal wave functions of N, A, C (color-octet particle) in the coordinate space, Xd 
describes the relative wave function of the deuteron, xaa and xcc represent the relative wave functions between 
As and Cs (in the single A A channel case, the CC component is absent), and C(S/) stands for the spin-isospin wave 
function of the corresponding system. It should be specially mentioned that in the form of such a wave function, 
normally called channel wave function [l3j], the totally anti-symmetric effect is implicitly included in the resultant 
relative wave function by solving the RGM equation and then projecting to the physical states. The channel wave 
functions of relevant systems are plotted in Fig. [2] 



































FIG. 2: Relatively wave functions in S— wave in the extended chiral SU(3) quark model: (a) for deuteron, (b) for xaa in the 
single A A channel case for d*( 2380), (c) for xaa in the coupled A A and CC channel case for d*( 2380). 

For the sake of convenience, we expand the relative wave function in the following: 

X( R ) = Xj Ci ex P 

We would also mention that the D -wave contribution is omitted due to its relevant smaller contribution, although 
both the S- and D -wave functions exist in our resultant wave functions. 

With these wave functions, two additional assumptions are employed in the estimation of the decay width in 
d* —> dim. One is associated to the energy denominators in Eq. ©, where the pole position is simply taken, as 
usually did in K-matrix approximation approach. The other one is related to the directions of the two outgoing pions. 
Since the experimental data show that the angle between two outgoing pions is almost zero, namely, the pions are 
propagating in the same direction, we can employ a bilinear condition for the momenta of these pions, k± || hj || z, so 
that the calculation can be much more simplified without loosing the major characters of such a process. 


III. NUMERICAL RESULTS AND DISCUSSION 

In the calculation, the masses of deuteron, A, nucleon and pion are taken from Particle Data Group 0- The mass 
of d* is Md * = 2Ma — e with e being 80 — 90 MeV for the single AA channel case and ~ 84 MeV for the coupled 
AA and CC channel case, respectively. The value of G is already fixed by using the A — > irN decay data. The decay 
width in the d* —> dim process with realistic wave functions from the RGM calculation can numerically be obtained 
by using Eq. ©• _ „ 

Moreover, the experimental data filial. I20l . l2ll| and one of theoretical calculations |22i] showed that for the d* resonance 
at y/s = 2.37 GeV, the decay cross section in the d* -A pmT + iT~ process is about 0.20 mb which is comparable with 
that of 0.24 mb in the d* —> dir 0 ir 0 process, and the decay cross section in the d* -A ppn°n~ process (also its mirrored 
channel d* -A nnir + ir°) has a visible value of about 0.10 mb as well. Therefore, contributions in these processes should 
also be accounted for in the d* width estimation. Using Breit-Wigner formulism and those cross section data, one 
estimated the branching ratios of various decay modes [1, 0]. For reference, we tabulate them in the second last 
column in Table [I] In order to consider the effect of isospin symmetry breaking of pions, namely the mass difference 
between and 7r°, we calculate the cross section in the d* -A dir + ir~ process explicitly. The obtained cross section 
in this calculation is about 1.83 times that in the d* -A dir°n a process, which is slightly larger than the multiple of 
1.6 estimated in Ref. 0,0- Based on resultant decay widths for d* -A dir°ir° and d* -A dir + iT and cross sections 
mentioned above, we get the branching ratios, and consequently the partial decay widths, for all possible decay modes. 
Finally, we achieve the total width of d*. 

In order to see the effect of the dynamical structure on the decay width, we calculate the width in the single AA 
channel case with e = 80 or 90 MeV (the corresponding mass of d* is about 2384 MeV or 2374 MeV, respectively) 
to compare with the result reported by Bashkanov et al. juj, where the decay width is reduced by the phase space 
effect only. We tabulate resultant decay widths for all the possible channels in Table Q] From this table, one sees that 
in the single A A channel calculation, no matter in which case (M^* = 2384 MeV or Md* = 2374 MeV), the resultant 
total width of d* justifies the fact that in a composite system, due to the binding behavior, namely the attractive 
interaction between constituents, the decay width of the system is much smaller than the total decay widths of its 
constituents if they were assumed to be free particles. And even more, deeper binding would cause narrower width. 
This feature is reasonable, because the width is not only related to the phase space, but also depends on the overlap 
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TABLE I: Decay width 



Ours 

Ex 

pt. 

21, 23] 

AA 

two channels AA + CC 

[4, 20,. 

M d , (MeV) 

2384 

2374 

2380 

2375 

channel 

T(MeV) 

T(MeV) 

B r 

T(MeV) 

B r 

T(MeV) 

d* —> d7T U 7T U 

22.6 

17.0 

13.3% 

9.2 

14(1)% 

10.2 

d* —» d7r + 7r“ 

41.5 

30.8 

24.3% 

16.8 

23(2)% 

16.7 

d* pmr u n u 

18.8 

14.2 

11.3% 

7.8 

12(2)% 

8.7 

d * —► pmv + n~ 

47.1 

35.4 

27.8% 

19.2 

30(4)% 

21.8 

d* —¥ pp7r u 7r” 

9.4 

7.1 

5.65% 

3.9 

6(1)% 

4.4 

d* —» nmv + n u 

9.4 

7.1 

5.65% 

3.9 

6(1)% 

4.4 

d* —> pn 

18.8 

14.2 

12.0% 

8.3 

12(3)% 

8.7 

Total 

167.6 

125.8 

99.9% 

69.1 

103(14)% 

74.9 


of the wave functions of the bound states d* and deuteron. In comparison with the estimated width of about 160 
MeV with the binding energy of 90 MeV by Bashkanov et al. 0, where the effect of the phase space is considered 
only, the contribution to the width from the dynamical structure of the system is about few tens of MeV. This tells us 
how important the effect of the dynamics on the width of an unstable composite system is, namely, the decay width 
is not only related to the phase space, but also depended on the dynamical structure of the system. It also shows that 
the width of d* in the single AA channel case where the mass of d* coincides the experimental data of 2384 MeV still 
far exceeds the experimental value of 75 MeV. This means that the AA structure alone cannot provide a reasonable 
width of d*. 

With the same scenario, we further exam the width contributed by the AA component in d* if d* has the AA+CC 
structure proposed in Refs. [M3- The results are also tabulated in Table |T] It shows that with the wave function 
of the AA component in Ref. m the decay widths for the d* —> d7r°7r° and d* — > dir + n modes are about 9 MeV 
and 17 MeV, respectively. If we further consider the d* — 1 prnnr, ppTr 0 TT~, nmr + ir a , and NN modes, the total width 
would be about 69.1 MeV. 

Here we would like to mention that in the RGM calculation, the trial wave function of the d* system is assumed 
to have two major components, AA and CC, which are totally anti-symmetrized. Solving the RGM equation, one 
obtains the relative wave functions of the system. By projecting the resultant wave function onto the cluster internal 
wave function in each component, we get the inter-cluster relative wave function, namely the channel wave function, 
for corresponding channel. Now, the contribution from the CC channel via the quark exchange is included in the 
projected wave function (or channel wave function) Xaa (1 ?) already [l3[. We should specially emphasis that the 
channel wave functions obtained in Eq. (7) are orthogonal to each other. Therefore, in the lowest order, by using this 
channel wave function, there is no quark exchange between the two physical particles and thus the colored clusters 
(color octet) cannot turn into the un-colored clusters (color singlet). As a consequence, the width contributed by the 
projected CC component would almost be zero. Combining this point with the contribution from the AA component, 
one sees that total width of d* in our A A +CC model is about 69.1 MeV, which is compatible with the experimental 
data of 75 MeV. Apparently, because the fraction of the wave function of the CC component in our AA+ CC model 
is about 67%, the resultant width of d* justifies our assertion that the d* resonance is a hexaquark-dominated exotic 
state. 

Finally, it should also be mentioned that the existence of d* should further be checked in other experimental 
processes. Now, except the q(or e)d reaction and pp collision, the strong decay of the hidden heavy flavor meson, such 
bb meson and cc meson, is also a place to hunt for d*. In particular, searching for its anti-particle d* in these processes 
is even plausible, because the anti-deuteron d and consequently d* can only be created from quark-pair productions, so 
that the background would be very clean (24|. Now, at yfs ~ 10.6GeV, the integrated luminosity is about 470/6” 1 at 
BaBar, and about 3/6“ 1 at CLEO. And both Collaborations have observed d production at yfs « 10.6 GeV |25i.l26j. 
Thus, one might search for d* in the T(nS) —> d* + p + n process. Moreover, the Belle Collaboration has collected 
even more data of about 1000 /6 _1 around yfs ss 10.6 GeV, and they might have the chance to observe the d and d* 
productions in the similar process. Also, BEPCII/BESIII has reached an integrated luminosity of 1/5“ 1 at yfs = 4.42 
GeV and 0.57/5” 1 at yfs = 4.6 GeV. It might be possible to detect d* in the e + +e~ —> d* +p + n process as well. If 
one could observe d* in the data set accumulated by BaBar, Belle, CLEO, and BEPCII/BESIII, it would definitely 
be helpful in confirming the existence of d*(2380) and its structure. 
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Appendix A: Spin-isospin part 

The spin matrix element in the calculation is 

T \ (~~lS ABmAB f~iS AB m AB ry^A m A (^^B m B 

s / J ^S A m A ,S B m B S' A m' A ,S' B m' B S A rn A ,l^S B m B ,lv 

= Y,^-) S ' ab ~ SaB+Sa+Sb ~ S ' a ~ S ' b ^'bSab , 3 2 2 3 i{ 5 ^ 

J SaB 1 J23 1 /-ijrnj /~iS' AB m' AB 

1 1 S'ab J J ' 

t 0 _ AiJ rn i (^S'ab^'ab 

where a = \/2 a + 1 . For the present process, the initial d* and final deuteron have Sab = 3 and S' AB = 1 , respectively. 
Moreover, A and nucleon have Sa = S B = 3/2 and S' A = S' B = 1/2, respectively. Therefore, in the present case 
J 23 = j = 2 is restricted. 

Moreover, one can deal with the isospin matrix element similarly. The isospins of A, nucleon, d* , and deuteron are 
3/2, 1/2, 0, and 0, respectively. Then j = 0 and j '23 = 1 are constraint for isospin part. 


(Al) 

S A 1 Sb S' b 1 

S' AB J23 j \ S A J23 Sab J 
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